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INTRODUCTION ; 

Tha Ontario Hiniatry of Environment in co-operation with 
tha Hiniatry of Solicitor General and Ontario Hydro, baa 
developed an emergency response program for tha release of 
radlo-nuclidaa In tha avant of an accident at any of the 
nuclear plants in Ontario. 

Two important pieces of Information, required for this 
program, are, tha wind field (i.e. wind speeds and 
directions downwind of the planta) and tha ground level 
concentration level fields. Generally measurements of winds 
are sparse and hence either extrapolation from a single 
aeasureaent point or meteorological forecaats are needed to 
generate the wind fields. This introduces a degree of 
subjectivity to the determination of tha direction in which 
the radlo-nuclidea propagate downwind after their release. 

It would be desirable to develop an objective scheme to 
simulate the wind fields in the vicinity of the nuclear 
power plants. We propose a three dimensional wind field 
model for this purpose. The model ie briefly described in 
the next section. Several criteria were used in the 
development of the modal. These are described below: 

a) Tha model has to include sufficient physics to 
account for the complex flow pattern* which are expected 
near a lakeshora. It is noted that all of Ontario's nuclear 
power plants are located on the lake shore. 

b) Tha model has to be computationally feasible to run 
in a short time to enable usto effectively respond in an 
emergency situation. 

The above criteria are very often In conflict with each 
other. The model described in this paper attempts to fulfill 
these criteria in an optimum manner. 
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The model Is tine dependent and is governed by the basic 
Navier Stokes equation. Atmosphere 1* assumed to be In a 
hydrostatic equilibrium I.e. vertical accelrations are 
neglected. Topography is included by defining a vertical 
coordinate which includes the local terrain height. Eddy 
exchange coefficients for momentum and heat are 
parameterized in the model ( for a detailed description of 
the model see Kaddukuri 1983). 

The model has ten Irregularly spaced vertical levels. 
Spacing is closer near the ground where higher resolution is 
required and la larger elsewhere. 

Figure 1 shows the model domain and Figure 2 the topography, 
as it la implemented for the Pickering Nuclear Generating 
Station. The model grid size la 5 km and covers an area 145 
km x 145 km. 

The model uses a preprocessor to compute the initial 
conditions that vary both horizontally and vertically. The 
preprocessor uses grldded (327 km) data (lsobaric heights 
and temperatures) provided by Canadian Heterological Centre 
(CHC) of Atmospheric Environment Services. Using orthognal 
polynomials the preprocessor interpolates this data down to 
5 km grid points. CMC balance winds or geostrophlc winds are 
used as initial winds in the model. 

Surface temperature which is the driving force for mesoscele 
circulations is derived from observed screen temperature, 
cloud cover and using empirical relationships developed by 
van Ulden and Hoistag (19B5). 

The preprocessor program writes all the Initial data onto 
appropriate files which are subsequently read by the main 
program. A time step of six minutes is used in the program. 

MODEL OOTFDTB IAH ETAHFLaH 

The model couuputes various parameters such as isobarlc 
heights, temperature, horizontal and vertical components of 
the wind, potential temperature etc. Data can be averaged 
over any desired time interval. Usually hourly averages are 
reported unless specified. Vertical averaging for wind speed 
is used to reduce noise in the data. Typically the model is 
used to forecast winds for six to twenty four period. 

Figure 3. shows the surface map for great lakes area on June 
24, 1988. Southern Ontario under the Influence of a high 
pressure system had clear skies and temeperatures were in 
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the uppper twenties. Figure 4. shows the initial conditions 
as computed at 8.00 AM local time. A« the surface 
temperature increases due to soair Insolation a lake breeze 



Figure 5. shows the model output at 15.00 PH. The 



tempera 

is set up. Figure 5. snows tne »»«■ mmnamw •> --■ 

shift In wind direction and magnitude la quite evident 

The mesoscale model Is also being verified with the data 
from a field study being conducted by Atomic Energy Board of 
Canada, Environment Canada, Ministry of Solicitor General 
and Ministry of Environment. One such comparison is 
presented in Figure 6. The agreement between observed and 
predicted wind speed and direction is quite good. 

nnHCLIIBIOHB I 

The model is capable of generating good mesoscale wind 
forecasts. Model verification with field data is currently 
underway. Initial results show a good agreement between 
observed and predicted windflleds. Additional tests are 
currently underway. After evaluation Is completed, the model 
will be implemented on microcomputer for use In emergency 
planning by Ministry of Solicitor General. 

s-utube FLAMS I 

After the model has been verified with the field data it 
will be implemented on a COMPAQ 386 microcomputer. Software 
will be developed to Interface directly with OME's HDAS 
(Heterological Data Acquisition System) so that 
meterological data can be transferred directly without 
manual entry. This total pacakage will be then be able to 
make mesoscale wlndflled forecasts on a routine basis and 
the data then could be used for emergency planning. 
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Fig. 1 i fteaoncala Hind Field Mattel Dcnain 
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Fig. 3 
SURFACE MAP 
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